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Abstract 

We compute spectral densities of momentum and R-charge correlators in thermal N = 4 Yang 
Mills at strong coupling using the AdS/CFT correspondence. For u ~ T and smaller, the spectral 
density differs markedly from perturbation theory; there is no kinetic theory peak. For large u>, 
the spectral density oscillates around the zero-temperature result with an exponentially decreasing 
amplitude. Contrast this with QCD where the spectral density of the current-current correlator 
approaches the zero temperature result like (T/u;) 4 . Despite these marked differences with pertur- 
bation theory, in Euclidean space-time the correlators differ by only ~ 10% from the free result. 
The implications for Lattice QCD measurements of transport are discussed. 
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I. INTRODUCTION 



The experimental relativistic heavy ion program has produced a variety of evidences which 
suggest that a Quark Gluon Plasma (QGP) has been formed at the Relativistic Heavy Ion 
Collider (RHIC) [1, 2]. The relative success of hydro dynamic approaches [3-7] suggests that 
the mean free path of the QGP is of order the thermal wave length ~ T [8]. This bold 
inference requires further theoretical and experimental corroboration. 

Theoretically, transport coefficients have been computed in the perturbative QGP using 
kinetic theory [9, 10]. However, it is important to obtain non-perturbative estimates for 
the transport properties of the medium. Kubo formulas relate hydrodynamic transport 
coefficients to the small frequency behavior of correlation functions of conserved currents in 
real time [11, 12]. In imaginary time correlation functions may be obtained from Lattice 
QCD measurements. Therefore, with a good model for the spectral density, one could hope to 
obtain a reasonable non-perturbative estimate for the transport properties of the QGP from 
Lattice QCD measurements. Recently, attempts to extract the shear viscosity [13, 14] and 
electric conductivity [15] have been made. More generally, the current-current correlator is 
being studied actively in order to extract interesting physical effects, e.g. dilepton emission, 
Landau damping, heavy quark diffusion, and in medium modifications of J/ip. [16-21]. 

Whenever there is a large separation between the transport time scales and the temper- 
ature T, Euclidean correlators are remarkably insensitive to transport coefficients [18, 22]. 
This has been studied in perturbation theory [22], where the transport time scale is l/(g 4 T), 
and in a theory with a heavy quark, where the diffusion time scale is M/T 2 or longer [18]. 
This insensitivity to transport is because the dominant transport contribution to the Eu- 
clidean correlator is an integral over the spectral density which is governed by a sum rule 
- the /-sum rule. Whenever, there is a quasi-particle description, the leading transport 
contribution to the Euclidean current-current correlator is proportional to the static sus- 
ceptibility times an average thermal velocity squared of the quasi-particle. Only with very 
accurate measurements and a solid understanding of the spectral density can more transport 
information be gleaned from lattice correlators. 

When there is no separation between the inverse temperature and transport time scales 
this conclusion requires further study. N = 4 Super Yang Mills (SYM) theory with a 
large number of colors and large 't Hooft coupling provides a theory where there is no scale 
separation and where the spectral density can be computed if the Maldacena conjecture 
is accepted [23]. This conjecture states that M = 4 SYM is dual to classical type IIB 
supergravity on an AdS x S$ background [23-25]. (See Ref. [26] for a review.) By solving the 
classical supergravity equations motion we may deduce the retarded correlators of strongly 
coupled TV = 4 SYM. 

Using the correspondence, Son, Policastro, and Starinets deduced the shear viscosity of 
strongly coupled M = 4 SYM [27]. Indeed the momentum diffusion constant is remarkably 
small, 7\ I (e+p) = 1/AttT. This result and subsequent calculations have led to the conjecture 
that the shear viscosity to entropy ratio is bounded from below, 7]/s > 1/Att [27-29]. The 
correspondence was also used to calculate the R-charge diffusion coefficient, D = 1/(2ttT) 
[30]. This transport coefficient is also remarkably small although there is no conjectured 
bound associated with this medium property. It is interesting to compute the full spectral 
densities of these correlators following the framework set up in these works. With the 
spectral density in hand, the feasibility of extracting transport from the Lattice in a strong 
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coupling regime can be assessed. 



II. STRESS TENSOR AND R-CHARGE CORRELATIONS 
A. Stress Tensor Correlations 

First let us recall the definitions [11]. The retarded correlators of T 00 and T° l are defined 
as follows 

Xee (k,u,) = dt /"d 3 xe +4 " x z^t)<[T 00 (x,t),T 00 (0,0)]> , (2.1) 

J —oo J 

/OO /* 
dt / d'xe+^'^ieit) ([T 0i (x,*),T°''(0,0)]) . (2.2) 
-oo J 

X g \ is then decomposed into longitudinal and transverse pieces 

Li uj / p uj \ 

X£(k, -) = ^ xS,(k, ") + [S- - — j x J 9 (k, a,) . (2.3) 

Through linear response, hydrodynamics makes a definite prediction for these correlators at 
small k and u; (see Appendix C for a review) 



^M-(e +P ) Jf^-^ t (2.4) 



Xj g (k,a;) = Vk \ k2 , (2.5) 
where e is the energy density, p is the pressure, rj is the shear viscosity, ( is the bulk viscosity, 

-ri+C 

and T s = -^j^- is the sound attenuation length. Using energy conservation we deduce 

uj 2 

Xee(k,w) = -j^X L gg- (2-6) 

Similarly we define the correlators of momentum fluxes 

/OO p 
dt d 3 xe +iujt - ik - x i9(t)([T il (x,t),T jm (0,0)}) . (2.7) 
-oo J 

Then using momentum conservation, we deduce that 

k i k l k i k m , u 2 



X 



X^(k,u,) = — xS.Ck.w), (2.8) 



- \ (> - *% m <* -) = Y 2 *- (k ' w) • (2 ' 9) 



The notation follows by assuming that k points in the x direction. 
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We will concentrate on the transverse piece of the momentum correlator. We define the 
spectral density as follows 

Imx^^k, uj) 



7T 



Using the functional form in Eq. (2.5) we deduce the Kubo formula 

7rp™(k,cj) 



lim lim 



V 



(2.10) 



(2.11) 



Euclidean correlators G yx . yx (pc, r 
density 

G v T xyx (k,r) = 



(T|, x (x, t)T^ x (0)) can be deduced from the spectral 
cosh (w (r - /3/2)) 



dw^(k,o;) 



(2.12) 



sinh (u/3/2) 

Further explanation of the Euclidean definitions and conventions adopted in this work may 
be found in [18]. 

The retarded correlator ^^(k, to) may computed using the AdS/CFT correspondence. 
The stress tensor couples to metric fluctuations h xy in the background of the black hole. 
We start with equation Eq. (6.6) of Policastro, Son and Starinets [30] for the fluctuation 
<j) = hy in the gravitational background. The linearized equations of motion for this field in 
the background are 



u 



to 2 



q 2 (l 



u 



uf 



up 



0. 



(2.13) 



Here <fik is the Fourier transform of with respect to four momentum k = (u, 0, 0, q) and 
we have defined to = u/(2nT), q = q/(2irT), and / = (1 — u 2 ). u — 1 corresponds to the 
the event horizon of the black hole; u = corresponds to spatial infinity. We will integrate 
this equation and set q = in what follows. 

Let us indicate the general strategy. Near the horizon u — 1, we develop the series 
solutions [31] 



u 



-in/2 



3w 2 _ -re , ■ tff s 

4 6 4 ^ 1 2 

1 + m 2 



(2.14) 

02 = 0T, (2.15) 

where <pi denotes the complex conjugate of 0i. The physical solution is <fii and not <f>\. 
Similarly near spatial infinity u = 0, we can develop a power series solution of this equation. 
These solutions are 



$ 2 



U 



-— log(u)$i + 







29 A 
-— m 


1 + W 2 u + 


G- 






72 J 



u 



(2.16) 
(2.17) 



We note that $i and $2 are real functions. Generally $j will be a linear combination X 
and 0i. Following [30], the physically correct solution of Eq. (2.13) is 0i which is a linear 
combination <3>! and $ 2 : 

1 = A$2 + -B$i. (2.18) 
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We may set A = 1 which is equivalent to the often assumed normalization condition 
4>i{u)\ u=Q = 1. Then the retarded Greens function is related by the correspondence to 

xT"(q,w) = lim/(<MA(0 fc )i , (2-19) 

where (4>k)i denotes the <f>i which solves Eq. (2.13) with k = (u;,q). Using the power series 
expansion near u « 0, and the fact that $i and $2 are real, we have the spectral density 1 

yxyxf \ Imy^^fu; O) 7T 2 iV 2 T 4 f 

Ptt { ] = Xtt { ,q) = ImB lim ^$ 2 <9 n ($i) , (2.20) 

vr 2 iV 2 T 4 

= 2 ImB . (2.21) 

4 V ; 

Thus the algorithm may divided into four steps: (1) Start at small u ~ and use the power 
series to find the initial condition and first derivative of $1 and $2- (2) Integrate $1 and 
$2 forward to determine the solution close to u ~ 1. (3) Using the power series close to 
u w 1 determine the linear combination of <p\ and <pi that is $1. Also do this for $2- (4) 
Then determine the linear combination of $1 and $2 that is <f>\. This fixes B in Eqs. (2.18) 
and (2.21), and consequently determines the spectral density at the frequency specified in 
Eq. (2.13). 

This algorithm was followed and the resulting spectral density is shown in Fig. 1(a). 
When the frequency tn is large we may follow this algorithm analytically and obtain the 
zero temperature result at large frequency. This is done in Appendix B and the resulting 
solution provides a good check of the numerical work. Further discussion of this figure in 
provided in Section III. 



B. R-Charge Diffusion 

A similar program may be followed to calculate the spectral density of the current-current 
correlator. First let us recall the definitions. Through linear response [11, 18], the diffusion 
equation predicts that the density- density correlator 

dt / d 3 xe + ^ ik - x z0(t) < [J°(x, t), J a °(0, 0)] ) (no a sum) , (2.22) 

•00 J 

has the following form at small k and o> 

y Dk 2 

"»-(*•") = =STW (2 ' 23) 

where Xs is the static R-charge susceptibility. The spatial current-current correlator is 
defined similarly 

/OO /' 
dt e + ^-^ie(t) ([Jl(x^), 4(0,0)]) (noasum), (2.24) 
-OO J 



1 As is common in the lattice community the spectral density is defined so that 
greens function differs by an overall sign from Rcf. [30] in order to conform with Ref. [18]. 



= 77. The retarded 

j=0 
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and may be broken up into longitudinal and transverse components 

X%<& ") = f - xL(k, u) + -g- X L rj (K w) • (2.25) 

The density-density correlator can be related to the longitudinal current- current correlator 

uj 2 k { W 

— XiVJV (k, u) = -J^XUK w) = x^. 7 (k, . (2.26) 

For k = there is no distinction between the longitudinal and transverse parts and therefore 
for k < T, x^(k,^) ^ xL(k,u;). 

The computation of these correlators at small k and u has already been performed by 
Policastro, Son, and Starinets [30]. From their computation (see Eq. 5.17b of that work) 
and the functional form of Eq. (2.23), the diffusion coefficient and static susceptibility are 

1 N 2 T 2 

D = , and Xs = • (2.27) 

2vrT' a 8 V ' 

To extend this computation to finite frequency it is simplest to compute Xjj(0> which is 
uncoupled from the other modes. Since k = 0, Xjj(0, w) is equal to Xjj(0> Perturbations 
in the R-charge current couple to the Maxwell field and the equations of motion for the 
Maxwell field in the gravitational background have been worked out. Consider the Maxwell 
field transverse to q, i.e. if q points in the z direction then = A x . The equation of 
motion for the Fourier components of in the gravitational background reads 

# + j& + ^(y-q 2 )& = 0. (2.28) 

with / = 1 — u 2 and as before k = (u,Q,Q,q), to = uj/(2hT) and q = q/(2,7rT). In what 
follows we set q = 0. 

The procedure mirrors the stress tensor computation. Near u = 1 we determine a power 
series solution 

01 = (1 -«)"** [! + ••■] , (2.29) 

02 = 0i. (2.30) 

0i is the physical solution and not 4>i- Near u = we determine the power series 

$i = wu + ... , (2.31) 
$2 = 1 + .... (2.32) 

Integrating from u « to u ~ 1 we determine the linear combination of $i and $2 that is 

01 

0! = $ 2 + B $! . (2.33) 

Then the correspondence states that 

N 2 T 2 

xL(0,o;) = __lim f(0_,OA(0 fc )i • (2-34) 
Since $2 and $1 are real and since $2 starts with one, the spectral density is 

irp^{0,u) _ ImXjAO,") N 2 T 



uj uj 16n 



Im B . (2.35) 
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FIG. 1: (a) The spectral density of the stress energy tensor, tt Ptt VX '(oj) / 'to normalized by the shear 
viscosity, r] AdS = ttN 2 T 3 /8. (b) The spectral density of the current-current correlator, irpjj/uj 
normalized by XsD = N 2 T 2 /16ttT. In both cases the dashed curves show the zero temperature 
results (Eq. (B16) and Eq. (B17)) normalized by the same factors. Due to a non-renormalization 
theorem in these channels, the zero temperature spectral densities in the free and interacting 
theories are equal [32, 33]. At finite temperature the kinetic theory peak does not exist in the 
strongly interacting theory. 



III. RESULTS 



The spectral density of the stress energy tensor 

p!™*(u) = dte + ^ J rf 3 x ([T«*(x,t), 7^(0,0)]) , (3.1) 

is shown in Fig. 1(a). Similarly the spectral density for the R-charge current- current corre- 
lator 

PjJ (Lo) = 7 ^J°°dte +i » t J d 3 x([j:(x,t),J^(0,0)]) (no a sum), (3.2) 
is shown in Fig. 1(b). These are normalized so that 

Kpjj{u) 



7] , and 



oj=0 



XsD, (3.3) 



ui=0 



where r] is the shear viscosity, Xs is the static R-charge susceptibility, and D is the R-charge 
diffusion coefficient. 

The remarkable feature of these spectral functions is the absence of any distinction be- 
tween the transport time scales and the continuum time scales. For comparison, consider the 
spectral density of the stress energy tensor in the free theory as worked out in Appendix A. 
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oa/(27iT) 



co/(2jiT) 



FIG. 2: The spectral density of the (a) stress energy tensor and (b) current-current correlators as 
in Fig. 1 but the zero temperature result has been subtracted and the absolute value taken. The 
plus or minus indicates the sign. The finite temperature spectral densities oscillate around the zero 
temperature result with exponentially decreasing amplitude. 

The full spectral density is a sum of the gauge, fermion, and scalar contributions 

(P^W = ( £ 7^) (l + 2n B M/2)) u\ (3.4) 

(pf/ X )f ermi on = (^) ^ + 7^2 (! " 2 r* M/2)) ^ » ( 3 - 5 ) 

(pf^Uiar = (4^) x ^ (a;) + 3^ (1 + 2^(^/2)) u\ (3.6) 

where for example (e + p) x is the partial enthalpy due to the free scalar fields - see Ap- 
pendix A. The free spectral function for the current-current correlator has a similar structure 
and is recorded in the appendix. Notice the delta function at the origin. According to the 
/ sum rule (see below), perturbations will smear the delta- function at the origin but will 
not change the integral under the peak. The peak indicates a large separation between the 
transport and temperature scales and is what makes kinetic theory possible. It is hopeless 
to apply a quasi-particle description to the strongly coupled Af = 4 case. 

Also shown in Fig. 1(a) and (b) is the spectral density at zero temperature for the stress 
tensor and R-charge correlators. At zero temperature, the strongly coupled spectral densities 
are identical to the free M = 4 spectral densities as required by a non-renormalization 
theorem [32, 33]. Fig. 2(a) and (b) show the absolute value of the difference between the 
finite temperature and zero temperature results on a log plot. The figure shows that the 
finite temperature result oscillates around the zero temperature result with exponentially 
decreasing amplitude. This is quite unusual and the exponential decrease is not due to any 
obvious Boltzmann factor. In QCD for instance, the spectral density of the current-current 
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FIG. 3: The Euclidean correlator for the (a) stress energy tensor correlator, Gtt VX , and the (b) 
current-current correlator, Gjj. The dashed curves show the free result for these Euclidean corre- 
lators. 



correlator approaches the zero temperature correlator as 2 , (T/u) 4 [35, 36]. The author has 
no explanation for the strong coupling results in M = 4 SYM. 

Finally, we may determine the corresponding Euclidean correlators by integrating the 
spectral density. For instance G^. yx is determined from p^. yx (u) 

G y T T{r) = J d 3 x(Tf(x,r)Tf(0,0)) , (3.7) 

sinhM/2) • (3 ' 8) 

The resulting Euclidean correlators for the stress tensor and current-current correlators are 
shown in Fig. 3(a) and (b) respectively. In each figure the strongly coupled case is compared 
to the free theory worked out in Appendix A. 

First consider the stress tensor correlations. Although the spectral density in the strongly 
coupled theory is markedly different from the free theory, in Euclidean space-time the cor- 
relators differ only by 10%. This 10% difference is significantly smaller than the errors 
associated with this correlator on the lattice [13, 14]. At least from the perspective of the 
AdS / CFT correspondence, it is hopeless to measure the transport properties of the medium 
in the T^ u channel where noisy gluonic operators dominate the signal. The figure also 
suggests that armed only with Euclidean measurements it is difficult to tell whether the 
theory is strongly interacting, i. e. whether there is a transport peak in the spectral density. 
Therefore, at least within the narrow framework of this work, rough agreement between 



2 There appears to be a discrepancy between two independent calculations; Ref. [34] finds that the difference 
falls as (T/uj) 2 . 
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quasi-particle calculations and lattice results [37, 38] says little about the validity of the 
quasi-particle picture. Only precise agreement can firmly establish validity. 

In weak coupling, the insensitivity to the transport time scale in the Euclidean correlator 
is readily understood from Eq. (3.8) and the free correlator Eqs. (3.4)-(3.6) [18, 22]. As 
perturbative interactions are turned on, the delta function at uo ~ is smeared by the 
inverse transport time scale. However, the integral under the peak remains the enthalpy 
times an average quasi-particle velocity squared, (v^/h). This is the physics content of the 
/-sum rule and is derived in Appendix D. When the smeared delta function is substituted 
into Eq. (3.8), the Euclidean correlator is roughly independent of the width of the peak, 
i.e. independent of the transport time scale. It seems that the strong coupling correlators 
remember this insensitivity to transport at weak coupling. 

For the current-current correlator the difference between the free and interacting cases 
is slightly larger, ~ 20%. Lattice measurements in the electromagnetic current-current 
channel are remarkably precise, ~ 0.5 - 2.0% [16, 17, 21]. This precision stems from the 
fact that the current-current channel correlates quark fields while the noisy T^ v channel 
correlates gluonic fields. With hard work and precise data, the AdS/CFT correspondence 
would suggest that some information about transport can be extracted from Euclidean 
measurements in this channel. However even in the current-current case, one needs a 
definite strategy to isolate the low frequency contributions. 

Note added. When this work was near completion, I learned that P. Kovtun and 
A. Starinets had nearly completed a similar computation [39]. 

Acknowledgments. I thank Peter Arnold, Raju Venugopalan, Andrei Starinets, 
and Peter Petreczky for useful discussions. D. Teaney was supported by grants from the 
U.S. Department of Energy, DE-FG02-88ER40388 and DE-FG03-97ER4014. 



APPENDIX A: THE FREE THEORY 



The free M = 4 Lagrangian is written as follows: 



C = 2tr | -^F 2 + h a (-ia ■ d)\ a - U.X^X^ , (Al) 

where "a" is a SU(4) index and "i" is a SO (6) index. Under flavor rotation, A a transforms in 
the fundamental representation of SU(4) and X 1 transforms as the fundamental representa- 
tion of SO(6). SU(4) and SO(6) are locally isomorphic. SU(4) matrices are parametrized as 
e ip A (T 4 ) A with trace normalization ti[(T A ) A (T 4 ) B ] = C A 5 AB and C 4 = 1/2. Similarly, SO(6) 
matrices are written as e %t3A<yT with trace normalization Cq = 1 [40]. The normalization 
convention adopted here has been fixed so that the AdS/CFT correspondence holds at the 
level of two point functions at zero temperature [32, 33]. 

Using the Noether method we compute the conserved R-charge current 

r A = 2 tr J {^j X i -iPx j + (T A ) ab \ a (?pj \ b 1 , (A2) 
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where d = d — d and gives symmetric Feynman rules. The spectral density pjj(u>) 
is easily computed at k = 0; the details of a similar computation appear explicitly in an 
appendix of Ref. [18]. The full spectral density is 



with 



P^H = (p^) fcrmion +(p^) scalar , (A3) 

/v°\ N 2 C 

(P^)scaiax = If J gg^r (l + 2n B M/2))o; 2 , (A4) 

/v°\ /V 2 *^ 

(^)fennian = Rf + ( l ~ ^ F {uj(3/2)) CU 2 . (A5) 



A 



Here (x° s ) x = \N 2 T 2 Cq and (x°) A = |7V 2 T 2 C4 are the free R-charge static susceptibilities 
associated with the scalars and fermions respectively. 

Similarly, we compute the spectral density of the stress tensor correlations. We first 
construct the canonical stress tensor using the Noether method and then construct the 
symmetric traceless Bellifante tensor as described in Weinberg's book [41]. The full stress 
tensor is written 

1 ~ V 1 Jgaugc ' V 1 Jfermion ' V 1 ^scalar ' l AU J 

with 

= 2 tr { F " pvK + f (-\ f2 ) } > ( A? ) 

= 2 tr + A + ^ Q* a H* • } ' ( A8 ) 

(^U = 2 tr {^A^X, + gT (-\d a X t d a x)j } . (A9) 

Now let us compute p y T x T yx . A straightforward (though lengthy) one loop computation 
calculation leads the spectral density p yxyx . The spectral density is a sum of the scalar, 
fermion and gauge boson contributions 



i%r% u) = (pf/*) gaugc + (pf r ^) fcrmion + (pf^) scalar . (mo) 

These contributions are 



(P V rr yx )„ = i 6 -^) <»8(u,) + ^- 2 (l + 2n B ( w 0/2)) ^ , (All) 
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e + 


V 
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e + 


V 



(P V rr yX ) {cm = ( -Y- ) + ^ (1 - 2 n F ( w /3/2)) W < , (A12) 



N 2 

'7 

(Pf/ X )scaia r = ) x uS M + ^2 (l + 2n B M/2)) (A13) 

where for example (e + p) x is the partial enthalpy due to the free scalar fields. Explicitly the 
partial enthalpies are (e +p) A = ^g-N 2 T 4 , (e + p) x = ±ff A 2 T 4 , and (e +p) x = ^-N 2 T\ 
Eqs. A7-A9 and Eqs. A11-A13 constitute the free spectral functions for the current- 
current and the tensor-tensor correlators. 
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APPENDIX B: WKB SOLUTION FOR LARGE to 



When tt> is large we may solve Eq. (2.13) by a WKB type approximation. We first 
introduce a change of variables 



\i l —^<t>{u) , (Bl) 



which obeys a Scrhodinger equation 

S + ^^)(- 3 + 4A + ^ + t ,')^0. (B2) 

This equation has two singular points at u = and u — 1. Away from the singular points 
we obtain the two WKB solutions, 

i>x ~ ^=cos(5(n) + 1 ), (B3) 
VP( U ) 

^2 ~ MS{u) + <h) , (B4) 

y/p[U) 

where p(u) and S(u) are the analogs of momentum and action and are given by 

P(u) = * 2 , , (B5) 



S(u) 



pU 

/ p(u') dv! — to tan _1 (-v/w) + tu tanlT^Vw) ■ (B6) 
Jo 



Following the WKB strategy we will solve the equation exactly near the singular points. 
Near u = we expand the potential and obtain the following equation 

**■( 3 .^,, = 0. (B7) 



du 2 \ Au 2 u 
The general solution to this equation is 

ij;(u) = Ci V / ttJ 2 (V4tt) 2 w) + C 2 VuN 2 (^/{4w 2 u) . (B8) 
Using the asymptotic expansions of the Bessel functions we obtain the matching formulas 

v^J 2 (2rav^) - 1= cos (S(u)-^) , (B9) 

>/av 2 (2ttV«) -> 1 -1 sm^)-^ . (BIO) 

In obtaining these formulas we have used the expansion of p(u) ~ Xo/y/u and S{u) ~ 2\x>y/u 
near zero. 

Near u = 1 we solve the equation 

S + 5<r 1 ^( 1+ra2 )^ - < B11 > 
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and obtain 

^ = C 1 (l-n)^ i f +C 2 (l-n)5+ l t . (B12) 

The physical solution (transformed according to Eq. (Bl)) contains only (1 — . Now 

it is a simple matter to show that 



u 



S(u) — — j + i sin (s(u) — 



Const (l-u)Hf. (B13) 



The constant is irrelevant to what follows. 

Putting together the connection formulas, Eqs. (B9), (BIO) and (B13), we find that 

uJ 2 (2tvVu) + iuN 2 (2tt}Vu) ^ Const (1 - u)"** . (B14) 

Up to a normalization constant this is the physical solution. The normalization is fixed by 
the requirement 4>{u)\ u=0 = 1 and the the series expansions, uN 2 (2tou) ~ ^ + 0(u) and 
uJ 2 {2xo^/u) f=s | ro 2 u 2 + 0(u 3 ) . Thus the physical solution in a neighborhood of u = is 

4>\{u) = -7rrD 2 nA r 2 (2rov / ^) + i7ito 2 uJ 2 {2tv^/u) . (B15) 
Substituting (f>i(u) into Eq. (2.19) we obtain the leading behavior 3 

^>)=6^ 4 " (BIG) 

This is the leading behavior at large frequency and agrees with the zero temperature result 
[24]. The WKB solution provides an excellent check of the numerical work. 

A similar WKB analysis applies to the R-charge correlator; the details are omitted. At 
large frequency the spectral density is 

N 2 

pJj{uj) = 32^^ ■ (B1?) 
APPENDIX C: HYDRODYNAMIC MODES 

In this appendix we the review the hydrodynamic modes to keep the treatment self 
contained. Consider a small perturbation from equilibrium. The stress tensor can be written 
as the equilibrium stress tensor plus small corrections 

<T 00 > = e + e(x,t), (CI) 
(T 01 ) = + ^(x,t). (C2) 

The velocity is small and is v = g/(e + p). 
The linearized hydrodynamic equations are 

d t e + d ig l = 0, (C3) 
dtgt + diT* = 0, (C4) 



In Fig. 1 this result is rewritten as irp(ui)/u> = ?/ AdS 7rtr> 3 . with r) AdS = (ttN 2 T 3 )/8. 



13 



where is 



T ij = 6 ij p _ v (gi v 3 + _ fidtv 1 ) - 5 ij (div l 



(C5) 



and we use metric (— , +, +, +) in this section. To solve these equations we first take spatial 
Fourier transforms, gk(t) = J c/ 3 xe~ ikx g(x, t) . Next we divide the momentum vector into 
transverse and longitudinal pieces, 



5k = gk + k #k 



where k • g£ = 0. For g£(t), the solution of these linearized equations is 

gkW = gk(0)e"^', 
where g£(0) is the initial condition. For e k (i) and g^(t) the solutions are 



g 2 1 aK 1 



Ir i.2 



e 2 



T.kH 



( ^^ + i^e k (0) 
c,e k (0) - «^k( C 



(C6) 

(C7) 
(C8) 



To connect these solutions with correlators we follow the framework of linear response 
[11, 18]. The definitions of the correlators used below are given in Section II A. We slowly 
turn on a small velocity field with a perturbing Hamiltonian 



H = H n 



tfW(x, t)T Uj (x,t) 



and switch it off at t — 0. i> l (x, t) obeys 

v*(x,t) = e et 9(-ty (x) . 
From the framework of linear response we have 

^<^(k,t)) = - x ^(k,tK(k). 



(C9) 



(CIO) 



(Cll) 



Writing v (k) = Vq (k) + kv^ (k), and substituting the tensor decomposition Eq. (2.3) into 
this equation, we obtain 



d t (g L (k,t)) = -x£(M)t#(k), 
$<g T (M)> = - X J 9 (k,t)v^(k). 

The velocity field can be eliminated in favor of the the initial values 

<g*(k,0)) = X2 sg (k)^ (k) , 

with the static susceptibility 

dte~ et / lxe- ik ' x <[T 0i (x,t),T 0j '(0,0)]) . 
11 



A s,gg 



(C12) 
(C13) 



(C14) 



(C15) 



As in Eq. (2.3) the static susceptibility is also broken up into longitudinal (Xs gg ) an d trans- 
verse {xj, gg ) components. 

Provided the wavelength of the perturbing velocity field is long, the system is in perfect 
equilibrium at time t = 0. The statistical operator describing this system at finite velocity 
is 



(C16) 



The stress energy tensor T^ v = (e + p) u^u u + pg^ u is 

<T 00 ) = e <T' y ) =p5 11 (T 0i ) =g l = ( e + p)v l , (C17) 
In this long wavelength limit then, the static susceptibility is from Eq. (C14) and Eq. (C17) 

dig') 



A s,gg 



{e + p) S 13 



(C18) 



or Xs,gg = Xs,gg = e + P- ft is a ^ so c l ear that the perturbing Hamiltonian does not change 
the average energy. This means that e(x, 0) = in Eq. (CI). 
From this discussion we first, eliminate v from the equations 



yy e _|_ p 

d t (g (k,0> = ~X flfl (M) e + p 



(C19) 

(C20) 



Comparing this result with the solution of the hydrodynamics equations given in Eq. (C6) 
and Eq. (C8) (with 6^(0) =0), we deduce the retarded correlators 



XggiKt) = (e + p)e 



[c s k sin (c s kt) + T s k 2 cos(c s A;t)] 



Xggfct) = rjk e 



2 -2£-t 

e+p 



To find the retarded correlator in frequency space, we integrate J °° e +tuJt and find 

iujT s k 2 — {c s k) 2 



Xgg(k,u) 
xj 9 (k,w) 



uj 2 — (c s k) 2 + iuT s k 2 
rjk 2 



-iuj + 



rjk 2 
e+p 



(C21) 
(C22) 

(C23) 
(C24) 



The spectral density is defined as the imaginary part of the retarded correlator by tt 

Pgg(^ W ) e + P ^sk 2 



UJ 



ix {uj 2 - c 2 s k 2 ) 2 + (ujT s k 2 ) 2 
1 T s k 2 /2 



e+p 



T s k 2 /2 



[n(uJ-c s k) 2 + (T s k 2 /2) 2 n(iu + c s k) 2 + (T s k 2 /2y 



PggiKu) 1 rik 2 

UJ 71 



UJ 2 + 



\e+Pj 



(C25) 
, (C26) 
(C27) 
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These results are one way to express the Kubo formulas. They can be expressed slightly 

differently using current conservation. Indeed taking imaginary parts of Eq. (2.8) and 
Eq. (2.9) we can express the Kubo formula in terms of the fluxes 

limlim^^ = ^ + C, (C28) 
fc^o uj 3 

Um Um np^(M (C29) 

uj^Ok-^O UJ 



APPENDIX D: THE / SUM RULE IN WEAKLY INTERACTING THEORIES 

For a weakly interacting theory the transport time scale Tr is much longer than the inverse 
temperature. The spectral density will have a sharp peak at small frequencies uj ~ 1/tr -C T 
which reflects this separation of time scales. We may derive a sum rule for this peak following 
the same steps and caveats as detailed at length in Ref. [18]. Fourier analysis together with 
Eq. (2.8) and Eq. (2.9) and yields the relations 

2T / A %~(k,^) « L ft;&(M)L 1/A , (Dl) 



UJ 



Here A is a cut-off that is much larger than the inverse transport time scale, but small 
compared to the temperature, l/r^ < A C T. The result of this integral does not depend 
A to first order in the scale separation. 

Since the time derivatives in Eqs. (Dl) and (D2) are evaluated at t ~ 1/A which is short 
compared to the collision time tr, the free streaming Boltzmann equation should be a good 
description of the initial equation of motion. In this effective theory 1/A may be taken to 
zero. 

To evaluate these time time derivatives of the retarded correlators, we consider the same 
perturbing Hamiltonian and setup as described in Appendix C - see Eqs. (C9)-(C13). At 
time t — the system is in perfect equilibrium with temperature T and velocity v (x). The 
thermal distribution function 

/0(X ' P) 5 e iz r -<iL>T Tl » U + Y f, (1 ± U) P 3 «<(*> • (M) 
with f p = 1/ (e Ep//T =F 1) . For short times the collision-less Boltzmann equation applies, 



d i d 

\-v — 

dt p dx i 



/(x,p,t) = 0. (D4) 



Here v p = p/E and must not be confused with the velocity field vq(x) and its Fourier 
transform t>o(k). The solution to this equation with the specified initial conditions is 

/(x,p,t) = / (x - v p t,p) . (D5) 
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Then the fluctuation in the momentum density is 

d 3 p 



(x,t) 



(2tt)* 



5f(x, p,t)p* 



(D6) 



with 5/(x, p, t) = /(x, p, t) — f p . Then taking spatial Fourier transforms with k conjugate 
to x and substituting the distribution function, Eq. (D5), we have 



d 3 P 



T J (2tt) 3 

For small times, we expand the exponential, and find 



-7p(i ± / p ) P y 



«g(k) . 



(D7) 



S*(M) 
with 

and 



e + p)^-it 2 A; 2 (e + p)(^) (5 



T(e + p) J (2tt 



d 3 p 

(2tt) 3 

d 3 p 



3 2,2 / \ /«p\ fc^" 



P- 



-2 „,2 



P" 



/p(l ± /p) ^ "f 



3 5 



Thus, from Eqs. (C12), (C13), (Dl), (D2), and (D8), we find 

2T £^p x T r x (k,co) » r(e+p)/| tt 

/•A j / ^2 

2T / — p^(k,o;) « T(e+p)(f 

. n 00 \ 5 



«S(k) 

(D8) 
(D9) 

(D10) 

(Dll) 
(D12) 
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